Available online at www.sciencedirect.com 
ScienceDirect 


Journal of Power Sources 163 (2007) 1080—1086 


JOURNAL OF 


www.elsevier.com /locate /jpowsour 


Short communication 


Simulation of abuse tolerance of lithium-ion battery packs 


Robert М. Spotnitz**, James Weaver*, Gowri Yeduvaka’*, 
D.H. Doughty °, Е.Р. Roth? 


* Battery Design LLC, 2277 DeLucchi Drive, Pleasanton, CA 94588, USA 
> Lithium Battery Department, Sandia National Laboratories, Albuquerque, NM 87185, USA 


Received 1 August 2006; received in revised form 1 October 2006; accepted 3 October 2006 
Available online 20 November 2006 


Abstract 


A simple approach for using accelerating rate calorimetry data to simulate the thermal abuse resistance of battery packs is described. The thermal 
abuse tolerance of battery packs is estimated based on the exothermic behavior of a single cell and an energy balance than accounts for radiative, 
conductive, and convective heat transfer modes of the pack. For the specific example of a notebook computer pack containing eight 18650-size 
cells, the effects of cell position, heat of reaction, and heat-transfer coefficient are explored. Thermal runaway of the pack is more likely to be 
induced by thermal runaway of a single cell when that cell is in good contact with other cells and is close to the pack wall. 


© 2006 Elsevier B.V. All rights reserved. 


Keywords: Lithium-ion battery; Thermal abuse; Battery pack; Accelerating rate calorimetry 


1. Introduction 


The safety of lithium-ion cells has been of foremost con- 
cern from their inception as evidenced by the number of safety 
devices developed specifically for lithium-ion cells [1,2]. How- 
ever, much less attention, at least in the technical literature, 
has been paid to the safety of lithium-ion packs. For example, 
lithium-ion packs for portable computers represent a significant 
portion of the market for 18650-size cells, and only a few papers 
(3,41 address the thermal behavior, much less the abuse toler- 
ance, of a typical pack. This short communication presents а new 
approach for estimating the thermal abuse tolerance of lithium- 
ion battery packs based on the behavior of individual cells. 


2. Model 


Approaches to modeling the thermal abuse of lithium-ion 
cells have been reviewed previously [5]. For pack modeling, the 
heat generated in the cell is determined experimentally and then 
used in an energy balance for the pack to predict the pack tem- 
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perature. The heat generation rate of the cell can be obtained by 
experimental characterization using accelerating rate calorime- 
try (ARC); ARC is widely used to characterize lithium-ion cells 
[6—8] but a review is beyond the scope of this short commu- 
nication. The ARC experiment provides direct measurement of 
the self-heating rate as a function of temperature. For simplic- 
ity, consider a hypothetical case where the ARC experiment for 
a single lithium-ion cell produces a normal curve (see Fig. 1). 
A more complex ARC profile produced from an actual exper- 
imental characterization of a cell could be considered, but an 
idealized situation is used here to help explain the method. The 
self-heating rate (SHR) is used to estimate the decomposition 
rate of the cell (Ка) by: 


SHR [K s7 !] 


= - 
ка В | Е Tmax [K] - Tmin [K] 


(1) 


where Tmin refers to the onset temperature for the thermal reac- 
tion and Tmax to the temperature at which the reaction is complete 
(see Fig. 1). The decomposition rate can be used to estimate the 
heat generated from the cell (qgen) by: 


gen [W] = На [Ла [871], 
Hixn = Meel [8]cp cen 08. КГ 1 (Tmax — Tmin) (2) 
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Fig. 1. Reaction rate as а function of temperature. 


Each cell is taken to be at a single temperature; that is, tem- 
perature gradients within cells are not considered. This is a 
reasonable approximation for small cells like the 18650 size 
[9]. Using this heat generation term, an energy balance can be 
written for each cell in the pack: 

сы! = 22 Boat + 3 + Qair,i + 4gen,i 


qwall,k,i = ое Aik raa(Ty — ТА) + hit Aik,contact(Tk — Т) 
comp, j,i = де Aijrad(T; = ТЁ) F hij Aij,contact(T; - Т) (3) 
Ча = X him Ага (Т) = Tj) 

m 


kik 


i. 
г Lik; F Lmkm 


where т; is the mass of the cell, ср, the heat capacity of the cell, 
o the Stefan—Boltzman constant, ¢ the emissivity of the cell, 
Aij raq the area of component j that is visible to cell i, hj the heat- 
transfer coefficient as defined above, k; the thermal conductivity 
of component i, and 1; is the distance from the center of mass 
to the surface of component i. A similar energy balance can be 
written for each wall that includes an additional term for heat- 
transfer from the wall to the external environment. The Battery 
Design Studio® software automatically computes the contact 
areas (Ajj contact) Of the various pack components as well as the 
line-of-sight areas (А; гаа) that are needed for the convective and 
radiative heat-transfer calculations, respectively [10]. The set of 
ordinary differential equations are solved using a Runge Kutta 
method. 

To apply these energy balances, a pack layout must be 
defined. For example, consider the pack layout shown in Fig. 2 
(based on the design presented by Maleki and Shamsuri [4]) 
with the thermal properties listed in Table 1. 


Table 1 
Thermal properties of pack 


Component Weight (g) Heat capacity Thermal conductivity 
dg! К!) (Wm! КТ!) 

Cell 43.8 0.67 n.a. 

Circuit board 50.0 0.7 500. 

Pack wall (x-z) 19.4 0.7 0.1 

Pack wall (х-у) 2.7 0.7 0.1 

Pack wall (у-2) 4.2 0.7 0.1 


insulated 


insulated 


Fig. 2. Pack design with eight cells. Three walls are insulated (y—x plane at z=0, 
x-z plane at y=0, and у-г plane at x=85 mm). 


The pack consists of eight cells, labeled as shown in Fig. 2. 
Three of the walls of the pack are posited to be insulated 
and treated as adiabatic: the bottom (х-2 plane), the right side 
(y-z plane), and the front (х-у plane). A small printed circuit 
board (PCB) is included in the pack, but spaced away from 
cell #3. 


3. Simulation results and discussion 


A common experimental method used to explore the abuse 
tolerance of a pack is to force one cell into thermal runaway 
and observe the resultant behavior: will other cells also go into 
thermal runaway or will the pack cool down? This can be an 
expensive test, so it is usually limited in scope. For example, 
only the cell in a single position might be brought into thermal 
runaway. However, with simulation, a wide number of exper- 
imental conditions can be explored. Here the effects of cell 
position (1—8), heat of reaction (corresponding to an adiabatic 
temperature rise of 270-330 К), and heat-transfer coefficient 
(10-100 W m~? K7!) are explored. 

To start the simulation, the cell in one position is set to a tem- 
perature that is high enough to induce thermal runaway (175 °C). 
Fig. 3 includes eight cases; each case shows the temperature pro- 
files of each cell in the pack as a function of time for the case 
where the adiabatic temperature rise of a single cell is 240 K, the 
initial pack temperature is 25 °С except for the one cell which 
is at 175 °С, and the pack heat-transfer coefficient to the envi- 
ronment is 10 W m~? K~!. In every case shown in Fig. 3, the 
pack does not go into thermal runaway, despite the low value 
of the heat-transfer coefficient; apparently the heat of reaction 
is too low to induce thermal runaway. If the pack is completely 
adiabatic, then all the cells go into thermal runaway and the pack 
reaches a steady-state temperature ~265 °C. 

Fig. 4 shows the case where the adiabatic temperature rise 
of a single cell is 300K and the heat-transfer coefficient is 
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Fig. 3. Temperature vs. time for each cell in the pack. Heat-transfer coefficient of 10 W m~? K7!. 
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Fig. 4. Temperature vs. time for each cell in the pack. Heat-transfer coefficient of 100 W m~? К-!. 
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Fig. 5. Temperature vs. time for each cell in the pack. Heat-transfer coefficient of 100 W m~? К-!. 
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Fig. 6. Temperature vs. time for each cell in the pack. Heat-transfer coefficient of 20 W m~? K~!. 
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100 W m~? K~!; in every case the pack goes into thermal run- 
away despite the large value of the heat-transfer coefficient. 

Fig. 5 shows the case where the adiabatic temperature rise of 
a single cell is 270 K and a high value of the heat-transfer coeffi- 
cient (100 W m~? K~!). In this case, runaway of cells 1—6 does 
not cause thermal runaway of the other cells in the pack. How- 
ever, if cell 8 goes into thermal runaway, cell 7 also goes into 
thermal runaway. Finally, if cell 7 goes into thermal runaway, 
then cell 8 goes into thermal runaway and, after about 10 min, all 
of the other cells. If the heat-transfer coefficient is reduced to a 
moderate value of 20 W m~? K~!, then the pack goes into ther- 
mal runaway no matter which cell goes into runaway first (see 
Fig. 6). This behavior can be rationalized as a balance between 
the ability of the pack to dissipate and absorb heat. 

The simulation results in Fig. 5 show that the pack is most 
susceptible to thermal runaway when cell 7 goes into thermal 
runaway. Runaway of cell 7 is most likely to induce thermal 
runaway of the pack because it has a large contact area with cells 
6 and 8. Cell 8 cannot dissipate much heat through the pack all 
and goes into thermal runaway. The additional energy provided 
by cell 8 is then able to drive cell 6 into thermal runaway. On 
the other hand, cell 6 is not as likely to drive the pack into 
thermal runaway because the heat from its exothermic reaction 
is distributed over more cells. 


4. Conclusion 


The simulation results presented here show that the thermal 
abuse tolerance of a pack is extremely sensitive to the exother- 


mic behavior of the cells. A small increase in the heat released 
by the exothermic reaction of a single cell can cause the pack 
to go into thermal runaway. This finding may explain why the 
number of safety incidents with laptop computers is increas- 
ing as the energy density of the cells increases. This work also 
shows the importance of heat-transfer from the pack and helps 
direct experimental work by identifying the cell most likely to 
induce thermal runaway in a pack. Most importantly, this work 
presents a new, simple approach for evaluating the abuse toler- 
ance of packs. 
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